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While time-dependent perturbation theory shows inefficient carrier-phonon scattering in semi- 
conductor quantum dots, we demonstrate that a quantum kinetic description of carrier-phonon 
interaction predicts fast carrier capture and relaxation. The considered processes do not fulfill en- 
ergy conservation in terms of free-carrier energies because polar coupling of localized quantum-dot 
states strongly modifies this picture. 
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I. INTRODUCTION 

Applications of semiconductor quantum dots (QDs) 
in optoelectronic devices rely on fast carrier scattering 
processes towards and between the discrete confined lev- 
els. These carrier transitions determine the dynamics 
of QD luminescence [ij or the operation of QD lasers 
[^|3{. For low carrier densities, where Coulomb scatter- 
ing can be neglected, carrier-phonon interaction provides 
the dominant scattering channel. In QDs only phonons 
with small momenta can efficiently couple to the confined 
carriers |j| . Then interaction with LA phonons does not 
contribute for large transition energies and only quasi- 
monochromatic LO phonons need to be considered. 

The simplest theoretical approach to electronic scat- 
tering processes is based on time-dependent perturba- 
tion theory. Fermi's golden rule for carrier transitions 
due to phonon emission or absorption contains a delta- 
function for strict energy conservation in terms of free- 
carrier energies of initial and final states and the phonon 
energy. When transition energies of localized QD states 
do not match the LO-phonon energy, efficient scatter- 
ing is inhibited (leading to the prediction of a phonon 
bottleneck) and only higher-order processes, like a com- 
bination of LO and LA phonons |j,|^, weakly contribute. 
Attempts to broaden the delta-function "by hand" imme- 
diately change the results |^ which underlines that this 
point should be addressed microscopically. The phonon 
bottleneck effect is still a debated topic, with experimen- 
tal evidence both for 0,11 and against it [TotlTlLll^. 

As in any coupled system, carrier-phonon interac- 
tion renormalizes both electronic and vibrational states. 
However, in bulk semiconductors or quantum wells with 
weak polar coupling, the net effect can be described by 
renormalized effective carrier masses, a small polaron 
shift of the band-edge, and lattice distortions only mod- 
ify the background dielectric constant for the Coulomb 
interaction of carriers. The broadening of the transition 
energies due to carrier-phonon interaction remains weak. 

For carriers in QDs, the discrete nature of localized 
electronic states changes the role of polaronic effects 
0, llSl [iJI ■ Restricting the analysis to a single QD state 



coupled to phonons, polaron effects can be obtained from 
an exact diagonalization of the Hamiltonian [1^ . While 
an extension to several discrete levels has been presented 
[l6j , the influence of the energetically nearby continuum 
of wetting-layer (WL) states, typical for self-assembled 
QDs, has not been included. Furthermore, only quasi- 
particle properties have been discussed which provide no 
direct information about the scattering efficiency for var- 
ious processes. Calculations of carrier transition rates 
based on the polaron picture are missing. 

We use a quantum kinetic treatment for carrier-phonon 
interaction in the polaron picture. As a first step, quasi- 
particle renormalizations due to the polar interaction for 
both QD and WL carriers are determined. For the QD 
states, the hybridization of one state with strong satel- 
lites of another state leads to a rich multi-peak structure. 
The WL states exhibit weak LO-phonon satellites. Cou- 
pling to the WL states provides a broadening mechanism 
for the QD states. 

Based on the spectral properties of QD and WL po- 
larons, quantum kinetic equations for the capture pro- 
cess (carrier transitions from the WL into the QD) and 
relaxation processes (transitions between QD states) are 
solved. For situations where, in terms of free-carrier ener- 
gies, energy conserving scattering processes are not pos- 
sible, the quantum-kinetic treatment provides efficient 
scattering rates. Even for the InCaAs material sys- 
tem with weak polar coupling, sub-picosecond scattering 
times are obtained. 



II. QUANTUM DOT POLARONS 

The single-particle properties of carriers under the 
influence of lattice distortions are determined by the 
retarded Green's function (GF), G^, which obeys the 
Dyson equation 
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Here a is an arbitrary (QD or WL) electronic state with 
free-carrier energy £„. In the polaron theory one usu- 
ally considers all possible virtual transitions from this 
state due to emission or absorption of phonons. This 
corresponds to a self-energy EJ^ for the carrier-phonon 
interaction where the population of the involved carrier 
states is neglected (electron vacuum). The correspond- 
ing retarded self-e nerg y in random-phase approximation 
(RPA) is given by 
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Assuming that the phonon system is in thermal equilib- 
rium, the phonon propagator (combined with the inter- 
action matrix elements) is given by 
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where monochromatic LO-phonons with the frequency 
ujlo are considered. The corresponding phonon popula- 
tion is given by n^o — l/{e'^'^^°/^'^ — 1) and the Frohlich 
interaction matrix element 
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contains the overlap between the electronic states and the 
phonon mode. For localized electronic states this acts as 
a form factor. The prefactor M|q = Aira-^iJlujLoY^'^ 
includes the polar coupling strength a and the reduced 
mass m. As a result of the above assumptions, the re- 
tarded GF itself depends only on the difference of time 
arguments and its Fourier transform can be directly re- 
lated to the quasi-particle properties. 

Due to energy separation between the discrete QD 
states and the WL continuum, polaronic effects in QDs 
are often computed by neglecting the presence of the WL 
0, Il3l [ij, [ig . For a single discrete level this amounts 
to the exactly solvable independent Boson model 15] 
and for several discrete levels it was shown to be nearly 
exactly solvable [I^l- In both cases, even for non-zero 
temperatures, the spectral function contains a series of 
sharp delta-like peaks. In real QDs, however, the in- 
teraction with the WL continuum (which might require 
multi-phonon processes) leads to a broadening of these 
peaks. The RPA accounts for this broadening effect while 
it retains a hybridization effect (see below) characteris- 
tic for the full solution. Therefore the RPA is expected 
to provide an adequate description in the presence of the 
continuum 19] . An additional source of broadening is the 
finite LO-phonon lifetime due to anharmonic interaction 
between phonons. 

For the numerical results presented in this paper we 
consider an InGaAs QD-WL system with weak polar cou- 
pling a = 0.06. The effective-mass approximation is as- 



sumed to be valid with rUe = 0.067mo for the conduc- 
tion band. For flat lens-shaped QDs the in-plane wave- 
functions of an isotropic two-dimensional harmonic po- 
tential are used while for the (strong) confinement in the 
direction perpendicular to the WL a finite-height poten- 
tial barrier is considered (see [la] for parameters and fur- 
ther details). To account for a finite height of the QD 
confinement potential, the calculations only include the 
(double degenerate) ground state and the (four-fold de- 
generate) first excited state, in the following called s and 
p-shell, respectively with s-p spacing and p-WL separa- 
tion of 40 meV. For the description of the WL states we 
use the following steps: (i) the WL states in the absence 
of QDs are described by plane waves for the in-plane part, 
multiplied by the state corresponding to the finite-height 
barrier confinement for the perpendicular direction, (ii) 
to describe the WL states in the presence of the QDs, 
the orthogonalized plane wave scheme, described in Ap- 
pendix A of Ref. [13], is used to construct WL states 
orthogonal to the QD states. Calculations are done for 
(3) a density of QDs on the WL n^ot = 10^° cm~^. Details 
on the calculation of the interaction matrix elements in 
Eq. Q with these wave functions for various combina- 
tions of QD and WL states can be found in Appendix 
B of Ref. jig]- Convergent results are obtained with 128 
points for the in-plane momentum radial integrals and 
50 points for the remaining angular integrations entering 
the interaction matrix elements. 

The Fourier transform of the spectral function, -2 Im 
G'^{uj), is shown in Fig. Q] for the fc = WL state and 
for the QD p- and s-shell (from top to bottom). In the 
absence of polar coupling to lattice distortions, the spec- 
tral functions are delta functions at the free-particle en- 
ergies indicated by the vertical lines. The dotted line in 
Fig. n^a) is the result for the fc = WL state without 
coupling to the QD states, indicating that their influence 
on the WL polarons is weak. The WL spectral function 
is broadened, the central peak exhibits a small polaron 
shift, and multiple sidebands due to LO-phonon emis- 
sion (absorption) appear to the right (left). The polaron 
broadening is a result of the irreversible decay in the con- 
tinuous WL density of states. 

The LO-phonon sidebands of the localized QD states in 
Fig.^b) and (c) are more pronounced and the hybridiza- 
tion of peaks from one shell with the energetically close 
sidebands of the other shell can be observed. This effect 
stems from the discrete nature of the localized states and 
requires that the coupling strength, which is modified by 
the form factors in Eq. I^J, exceeds the polaron damp- 
ing. If only coupling matrix elements diagonal in the 
state index would be considered, a series of sidebands of a 
state with discrete energy, spaced by the LO-phonon en- 
ergy, would be obtained. Off-diagonal coupling elements 
alone would lead to a hybridization of discrete levels as, 
e.g., in quantum optics where instead of the phonon- field 
a monochromatic light field coupled to a two-level sys- 
tem is considered 4]. When the level splitting equals 
the LO-phonon energy, in the limit of weak damping the 




FIG. 1; (a) Spectral function of electrons for the lowest WL 
state at fe = under the influence of polar coupling in the 
combined QD-WL system (solid line) and without coupling 
to the QD states (dotted line). For electrons in the QD p- 
shell (b) and s-shell (c) full coupling between all states (solid 
line) is compared to the case without coupling to other QD 
states (dotted line). Vertical lines show the corresponding 
free-carrier energies. Energies are given relative to the con- 
tinuum edge Eg in units of the phonon energy Tujjlo- The 
temperature is 300 K. 



splitting of each line is determined by the carrier-phonon 
coupling strength. Due to the interplay of diagonal and 
off-diagonal interaction matrix elements, the QD spectral 
functions in Fig. ^^b) and (c) show a series of satellites, 
each of them reflecting the hybridization. The asymme- 
try stems from the difference between level spacing (40 
meV) and LO-phonon energy (36 meV). The broaden- 
ing of peaks stems mainly from the coupling to the WL 
states. A finite LO-phonon lifetime of 5 ps due to anhar- 
monic interaction between phonons has been included in 
the calculations. 



III. CARRIER KINETICS OF RELAXATION 
AND CAPTURE PROCESSES 

In this section we study consequences of the renormal- 
ized quasi-particle properties on the scattering processes. 
Fermi's golden rule, which has been frequently used in the 



past, describes only transition rates from fully populated 
initial into empty final states. Proper balancing between 
in- and out-scattering events, weighted with the popula- 
tion / of the initial states and the blocking 1 — / of the 
final states, leads to the kinetic equation 
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Quasiparticle and Markov approximation can also be ap- 
plied to renormalized polaronic states which results in 
a rate-equation description for the population of these 
states 20] . 

A quantum-kinetic approach extends Eq. (jSJ in the 
sense that the delta-functions with free-carrier energies 
are replaced by time-integrals over polaronic retarded 
GFs. Furthermore, the population factors are no longer 
instantaneous but explicitly depend on the time evo- 
lution. This is the time-domain picture for the inclu- 
sion of renormalized quasi-particle properties (beyond 
a quasi-particle approximation and beyond Markov ap- 
proximation). Using the generalized Kadanoff-Baym 
ansatz (GKBA) [2ll, the quantum- kinetic equation has 
the form 
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The phonon propagator D^ follows from Eq. ^ by re- 
placing T — > — T. A Markov approximation in the renor- 
malized quasi-particle picture corresponds to the assump- 
tion of a slow time-dependence of the population fa{t2) 
in comparison to the retarded GFs such that the popu- 
lation can be taken at the external time ti. The Boltz- 
mann scattering integral of Eq. (O follows if one addi- 
tionally neglects quasi-particle renormalizations and uses 
free-carrier retarded GFs |22| . 

To demonstrate the influence of quantum-kinetic ef- 
fects due to QD-polarons, we first study the relaxation 
of carriers from p-shell to s-shell for the above discussed 
situation where the level spacing does not match the LO- 
phonon energy such that both, Fermi's golden rule and 
the kinetic equation (O predict the absence of transi- 
tions. A direct time-domain calculation of the polaron 
GFs from Eq. Q-Q together with Eq. © is used. We 
assume an initial population fsito) — 0, /p(io) = 0.3 and 
start the calculation at time to. While this example ad- 
dresses the relaxation process itself, more advanced cal- 
culations would also include the carrier generation via op- 
tical excitation or carrier capture discussed below. Then 
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FIG. 2; Temporal evolution of the QD population due to 
carrier-phonon scattering between p-shell and s-shell for an 
energy spacing larger than the LO-phonon energy. The solid 
lines correspond to a quantum-kinetic calculation whereas for 
the dotted line the Markov approximation is used together 
with polaronic spectral functions. 
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FIG. 3: Time- evolution of the QD p-shell (solid line) and s- 
shell (dashed line) electron population due to carrier capture 
from the WL including the effect of carrier relaxation between 
QD shells. If only direct capture processes are considered, the 
dashed-dotted and dotted lines are obtained for p-shell and 
s-shell, respectively. 



ambiguities due to initial conditions can be avoided since 
the population vanishes prior to the pump process which 
naturally provides the lower limit of the time integral 
in Eq. jnj. In practice, we find that within the GKBA 
results weakly depend on the details of the initial condi- 
tions. 

The evaluation of the quantum-kinetic theory (solid 
lines in Fig. |21) yields a fast population increase of the 
initially empty QD s-shell accompanied by oscillations 
which reflect in the time-domain the hybridization of cou- 
pled carrier and phonon states. The analogy to Rabi os- 
cillations has been pointed out in Ref. [3|. If one uses 
the Markov approximation together with polaronic re- 
tarded GFs in Eq. ||HJ), such that quasi-particle renor- 
malizations are still included, these transient oscillations 
disappear. In both cases the same steady-state solu- 
tion is obtained which corresponds to a thermal pop- 
ulation at the renormalized energies. The equilibrium 
solution can be obtained from the polaron spectral func- 
tion using the Kubo-Martin-Schwinger (KMS) relation, 
/a = -/^ /l'^) ImG^(w) where /(w) is a Fermi 
function with the lattice temperature. Note that par- 
ticle number conservation is obeyed in Fig. |21 since the 
degeneracy of the p-shell is twice that of the s-shell. 

Another important process is the capture of carriers 
from the delocalized WL states into the localized QD 
states. For the used QD parameters, where the spacing 
between the p-shell and the lowest WL state (40 meV) 
exceeds the LO-Phonon energy, again Fermi's golden rule 
and Eq. (jSJ predict the absence of electronic transitions. 
For the numerical solution of Eq. @ we use now as initial 
condition empty QD states and a thermal population of 
carriers in the polaronic WL states (obtained from the 
KMS relation) corresponding to a carrier density 10^^ 
cm~^ and temperature 300 K [23 ■ The dashed-dotted 
and dotted lines in Fig. O show the increase of the p- 



and s-shell population, respectively, when only capture 
processes are considered (scattering from a WL-polaron 
to a QD-polaron state due to emission of LO-phonons). 
Also in this situation the quantum-kinetic theory predicts 
a fast population of the initially empty p-shell. Albeit 
the large detuning (exceeding two LO-phonon energies) 
the direct capture to the s-shell is still possible but con- 
siderably slower. When both, direct capture of carriers 
as well as relaxation of carriers between the QD states 
are included in the calculation, the solid (dashed) line 
is obtained for the p-shell (s-shell) population. While 
faster capture to the p-shell states leads at early times 
to a p-shell population exceeding the s-shell populations 
(see inset of Fig.Ol, the subsequent relaxation efficiently 
populates the s-states. Since the WL states form a quasi- 
continuum, beating at early times is strongly suppressed. 
With the results in Figs. lll^l we have demonstrated the 
ultrafast (subpicosecond) carrier relaxation and fast (pi- 
cosecond) carrier capture for a material with weak polar 
coupling and 10% detuning between the transition ener- 
gies and the LO-phonon energy. This detuning is on the 
one hand sufficiently large for the alternative LO-I-LA 
mechanism proposed by Inoshita and Sakaki Q to fail 
and on the other hand small enough to illustrate the hy- 
bridization of one state with sidebands of the other states. 
We find that the fast scattering is not related to the near 
resonance condition and in fact relatively insensitive to 
the detuning between transition energies and LO-phonon 
energy. The spectral functions of the coupled QD-WL 
system for various detunings, ranging from resonance to 
a large mismatch of 40%, are shown in Fig. ^ For bet- 
ter visibility only the curves for the s- and p-shells are 
displayed, while the WL-states are also included in the 
calculations. As seen in Fig. ^ the spectral function for 
the WL states is only weakly infiuenced by the coupling 
to the QDs. In all three cases of Fig. ^ there is sub- 
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FIG. 4: (a) Spectral function of electrons in the coupled QD- 
WL system for the p-shell (solid line) and s-shell (dashed line) 
using various energy spacings A_B between s-shell, p-shell and 
WL {k = 0) in units of the LO-phonon energy Kojlo- Vertical 
lines indicate the positions of the unrenormalized QD states, 
Eg is the continuum edge of the WL states. The temperature 
is 300 K. 



stantial overlap between the s-shell and p-shell density of 
states which points to efficient transition processes. This 
overlap is due to the multi-peak-structure which contains 
the series of phonon sidebands spaced by the LO-phonon 
energy and their hybridization. From top to bottom in 
Fig. ^ the peak splitting increases with detuning. 

The corresponding results for the carrier relaxation, as 
in Fig. [5] but for different detunings A_E, are shown in 
Fig. [S] The fast carrier relaxation towards an equilib- 
rium situation is retained in all three cases. The main 
difference is in the oscillation period, which is reduced for 
larger detuning due to the increased splitting in Fig. ^ 

A stronger influence of the detuning between transition 
energies and the LO-phonon energy is found for the cap- 
ture of carriers from the WL into the QD states. As can 
be seen in Fig. from the resonance situation to a de- 
tuning of 40% the capture efficiency is reduced by about 
one order of magnitude. Nevertheless, a significant occu- 
pancy can be reached within several ten picoseconds. The 
reduced capture efficiency is related to a reduced overlap 
between the WL and QD spectral functions for increas- 
ing detuning (which is mainly because the WL states are 
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FIG. 5: Temporal evolution of the QD population due 
to carrier-phonon scattering between p-shell (initially pop- 
ulated) and s-shell (initially empty) for different energy spac- 
ings A_B corresponding to Fig. 0] 
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FIG. 6: Time- evolution of the QD p-shell (top) and s-shell 
(bottom) electron population due to carrier capture from 
the WL including the effect of carrier relaxation between 
QD shells for different energy spacings AiJ corresponding to 
Fig.Bl 



weakly influenced by the QD states). In comparison to 



this, the strong interaction between s- and p-states main- 
tains a strong overlap between their spectral functions. 
As a consequence the relaxation time is less sensitive to 
the detuning. 

In summary, the quantum-kinetic treatment of carrier- 
phonon interaction explains the absence of a phonon 
bottleneck in terms of scattering between renormalized 
quasi-particle states. A quasi-equilibrium situation is 
reached on a ps-timescale at elevated temperatures even 



in materials with weak polar coupling. 



Acknowledgments 

This work was supported by the Deutsche Forschungs- 
gemeinschaft and with a grant for CPU time at the NIC, 
Forschungszentrum Jiilich. 



[1] D. Morris, N. Perret, and S. Fafard, Appl. Phys. Lett. [16] 

75, 3593 (1999). 
[2] P. Bhattacharya, K.K. Kamath, J. Singh, D. Klotzkin, [17] 

J. Phillips, H.-T. Jiang, N. Chervela, T.B. Norris, T. 

Sosnowski, J. Laskar, M.R. Murty, IEEE Transactions 

on Electron Devices 46, 871 (1999). 
[3] G. Park, O.B. Shchekin, D.L. Huffaker, and D.G. Deppe, [18] 

Electron. Lett. 36, 1283 (2000). 
[4] T. Inoshita and H. Sakaki, Phys. Rev. B 56, 4355 (1997). [19] 

[5] T. Inoshita and H. Sakaki, Phys. Rev. B 46, 7260 (1992). 
1^6] H. Jiang and J. Singh, IEEE J. Quantum Electron. 34, 

1188 (1998). 
[7] J. Urayama, T. B. Norris, J. Singh, and P. Bhattacharya, 

Phys. Rev. Lett. 86, 4930 (2001). [20] 

[8] A.W.E. Minnaert, A.Yu. Silov, W. van der Vleuten, 

J.E.M. Haverkort, and J.H. Wolter, Phys. Rev. B 63, [21] 

075303 (2001). 
[9] S. Xu, A.A. Mikhailovsky, J.A. Hollingsworth, and V.I. 

Klimov, Phys. Rev. B 65, 045319 (2002). [22] 

[10] E. Tsitsishvili, R. v. Baltz, and H. Kalt, Phys. Rev. B 

66, 161405 (2002). 
[11] E. Peronne, F. Fossard, F.H. Julien, J. Brault, M. 

Gendry, B. Salem, G. Bremond, and A. Alexandrou, 

Phys. Rev. B 67, 205329 (2003). 
[12] F. Quochi, M. Dinu, L.N. Pfeiffer, K.W. West, C. 

Kerbage, R.S. Windeler, and B.J. Eggleton, Phys. Rev. 

B 67, 235323 (2003). [23] 

[13] K. Krai and Z. Khas, Phys. Rev. B 57, 2061 (1998). 
[14] O. Verzelen, R. Ferreira, G. Bastard, T. Inoshita, H. 

Sakaki, phys. stat. sol. (a) 190, 213 (2002). 
[15] G.D. Mahan, Many- Particle Physics (Plenum, New York, 

1990). 



T. Stauber, R. Zimmermann, and H. Castella, Phys. Rev. 
B 62, 7336 (2000). 

In general, Green's functions and self-energies depend on 
two quantum numbers a, (5. For the situation described 
in this paper, off-diagonal contributions can be neglected, 
which has been verified numerically; see also |j, [13 . 
T.R. Nielsen, P. Gartner, and F. Jahnke, Phys. Rev. B 
69, 235314 (2004). 

Comparison with an alternative method supports this 
statement: the first term of a cumulant expansion, which 
is known to be exact for the independent Boson model 
pa j. agrees well with the RPA for quasi-continuous elec- 
tronic states. 

O. Verzelen, R. Ferreira, and G. Bastard, Phys. Rev. B 
62, 4809 (2000). 

for a review see: H. Haug and P. Jauho, Quantum Kinet- 
ics in Transport & Optics of Semiconductors, Springer- 
Verlag, Berlin, 1996. 

The quantum-kinetic equation provides energy- 
conserving scattering processes only in the long-time 
limit. The corresponding spectral function in Fig. Q is 
obtained from a complete Fourier transform whereas 
on a short timescale the peaks shown there are only 
partially emerged (since only a finite-time contribution 
of the retarded GF enters in the quantum-kinetic 
equation at a given time.) 

For small carrier densities, we use the polaron GF for 
the electron vacuum since it is weakly influenced by pop- 
ulation effects, see P. Gartner, L. Banyai, and H. Haug, 
Phys. Rev. B 60, 14234 (1999). 



